APPROACHES TO THE MONOPOLE-DYNAMIC DIPOLE 
VACUUM SOLUTION CONCERNING THE STRUCTURE OF 
ITS ERNST'S POTENTIAL ON THE SYMMETRY AXIS 
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Abstract The FHP algorithm [1] allows to obtain the relativistic multi- 
pole moments of a vacuum stationary axisymmetric solution in terms of coefficients 
which appear in the expansion of its Ernst's potential £ on the symmetry axis. 
First of all, we will use this result in order to determine, at a certain approxima- 
fNj 1 tion degree, the Ernst's potential on the symmetry axis of the metric whose only 

multipole moments are mass and angular momentum. 

By using Sibgatullin's method [2] we then analyse a series of exact solu- 
te. ' tions with the afore mentioned multipole characteristic; besides, we present an 
O ' 

dimensionless parameter. The calculation of its multipole moments allows us to 
Q\ . understand the existing differences between both approximations to the proposed 

pure multipole solution. 
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1. INTRODUCTION 



As it is well known, the relativistic multipole moments denned by Geroch [3] 
and Hansen [4] allow us to characterize, in an unique manner, vacuum stationary 
metrics. Particularly, Schwarzschild's solution (spherically symmetric) can be 
described as that one whose unique multipole moment is monopole. 

In [5] we have introduced a static and axisymmetric solution of the Ein- 
stein vacuum field equations whose only multipole moments are mass and the 
quadrupole moment. Our aim in the present work is to generalize the above 
mentioned result to the non-static case, by searching the stationary and axisym- 
metric vacuum solution which in addition to mass possesses only a dynamic dipole 
moment. 

Fodor, Hoenselaers and Perjes [1] have developed an algorithm to calculate 
multipole moments in terms of the coefficients arising in the expansion of the 
Ernst's potential £ on the symmetry axis as a power series on the inverse of the 
Weyl's coordinate z. According to this result it is possible to determine such coef- 
ficients corresponding to a solution with the mentioned multipole characteristics. 

There is a method, due to Sibgatullin [2], which allows to generate exact 
solutions of vacuum, stationary and axisymmetric field equations from the Ernst's 
potential on the symmetry axis, whenever its structure is a polynomial ratio. This 
method has been broadly used [6], [7], [8] and it has been recently completed by 
the introduction of some general expressions which make its use much easier [9]. 

Our first aim, is to prove that a finite set of coefficients which appear in the 
expansion of the potential £ on the symmetry axis is sufficient enough to describe 
a potential E = (1— £)/(l + £) of rational type on the symmetry axis. Therefore, 
the use of Sibgatullin's method, makes it possible to construct an Ernst's potential 
E in terms of these coefficients. Nevertheless, if we look for an exact solution of 
the M — J type, then the result of this condition leads to the Ernst's potential on 
the symmetry axis by means of a ratio of series. Hence, there is no finite number 
of coefficients which describe the Ernst's potential of Monopole-Dynamic dipole 
solution as a polynomial ratio. 

We then approach the M — J solution as the limit of a sequence of exact solu- 
tions which possess a progressively smaller number of multipole moments higher 
than the dipole. In spite of this, the behaviour of the multipole moments shows 
that the mentioned series approaches the M — J solution in a rather unexpected 
way. 

In order to obtain an alternative approach to the M — J solution, we will 
propose a series of approximate solutions described by the partial sums of the 
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expansions of the Ernst's potential £ as a power series in a dimensionless parameter 
J. To do that, we use Schwarzschild's solution as a seed solution and solve Ernst's 
equation at successive orders in the parameter and impose the corresponding 
structure of the M — J solution on the symmetry axis. 

As we will see, the order of magnitude of its multipole moments decreases as 
the order of approximation rises. That leads us to conclude that the parameter 
J controls the successive solutions and contributes to give a physical meaning to 
the approximate solution. 

A more realistic stationary solution would be a M — Q — J solution, i.e., the 
one having only mass, massive quadrupole moment and dynamic dipole, since a 
rotating object flattens and, hence, all its massive multipole moments represent 
such deviation from sphericity. Nevertheless, we can imagine so rigid an object 
that the M — J solution itself would be physically relevant. Besides, since the 
static case is already solved by M — Q solution, we want to discuss M — J solution 
and consider the M — Q — J solution as a generalization of both solutions. 



2. STRUCTURE OF THE M-J SOLUTION ON THE SYMMETRY 
AXIS 

Let us be £ the Ernst's potential of a stationary axisymmetric solution of 
vacuum Einstein's field equations [10] 

(er-i)A£ = 2£*(vo 2 , (i) 

1 - E 

being £ = — , where the Ernst's potential E is the complex function whose 

1 + E 

real part represents the norm of the Killing vector describing stationarity. On 
symmetry axis, this potential £ can be expanded by means of a power series of 
the inverse Weyl's coordinate z as follows: 

oo 

ap = 0,z) = Y,m n z- (n+1 ^ (2) 

n=0 

where p represents the Weyl's radial coordinate. 

Fodor, Hoenselaers and Perjes [1] have developed an algorithm which allows 
to calculate the Geroch [3] and Hansen [4] relativistic multipole moments, related 
to a vacuum stationary axisymmetric solution, in terms of the coefficients m n 
arising in the previous expansion (2). Both the result obtained up to multipole 
order 10 by the afore mentioned authors, and the calculations we have carried 
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out up to order 20 lead us to show that the relation between multipole moments 
and coefficients m n is triangular. That is to say, the multipole moment and the 
corresponding coefficient m n at every order, differ in a certain combination of 
lesser order rrik coefficients. Therefore, these relations enable us to determine 
unequivocally the coefficients rrih which allows to outline the expansion of £ in 
terms of the known multipole moments for any given solution. 

So we have obtained that the solution having only massive monopole and 
dynamic dipole is characterized by an Ernst's potential £ whose expansion on the 
symmetry axis provides the following coefficients m n up to the order 20, 

mo = M , mi = iJ 
mi = , ms = 

1 9 

m 4 = -MJ 2 

1 • 73 

m 5 = -— U 

m 6 = — M 3 J 2 
b 21 

m 7 = -i|M^ 

- = ll M5j2 ~ ^ MJ4 

m9 = -^M*J 3 i-—J 5 i 
9 3003 3003 

mi0 = Am 7 j 2 - — m 3 j 4 

10 429 7007 

39 357357 

mi2 = -J— M 9 J 2 - -i^-M 5 J 4 — — MJ 6 

143 119119 2909907 

43 ,rfr3. 265 „4r5. 13051 T 7- 

Wl3 = -2431 M J ? " 108927 M J i " 2(3369349 J z 

m 14 = M 11 J 2 - -i^M^ 4 + 1129 M3J 6 
221 20369349 10968111 

53 „, 10t3 . 10954 831513 „ ~ 

mis = M 10 J 3 i M 6 J 5 z M J'i 

5 4199 2263261 364385021 
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mi 6 

mi 7 



1 ,.13 T 2 40346 lrfl 4 4 5 4 „ 5t6 2419504 „ » 
-M 13 J 2 — — M 9 J 4 ■ M 5 J 6 -MJ 8 



323 6789783 15954939 16397325945 

3 ^12 x3- 12480070 ~. 293822614 „ . 
-M 12 J 3 i M 8 J°i M 4 J 7 i- 



323 1717815099 60123528465 

35634548 t9 



147575933505 



mis 



5 „i 5t2 862123 „ , 1t4 10703470 - t6 
-M 15 J 2 — -M J — — -— M 7 J - 



2261 245402157 36074117079 

11186023022 „ tS 



21103358491215 



173594465986 ,, fir7 99041865574 T „ ,, 2 365 ,,14 ^ 
= - 211033584 9121 5 M J ' " 87428199463605 " 52003 M J " 

6808829 _ Ml0j5 . 



736206471 



"7 2 = 



5 ^17 x2 6540151 , 3t4 9908406983 „ ~ t6 
-M 17 J 2 M 13 J 4 ■ M 9 J - 



3059 3681032355 21103358491215 

4674899812546 - t8 985078066594 „ , n 
-M 5 J 8 MJ 10 



4283981773716645 18971919283602285 

(3) 

where M and J represent Mass and Angular Momentum respectively. 

These expressions suggest that the coefficients m n of the potential £ repre- 
senting the Monopole-Dynamic dipole solution can be written in the following 

way 

(2fe+l)/4 

m 2k = M 2k+1 J 2 ™ °( 2n > 2k ) 

n=l 

, (fc = 0,l,...) (4) 

fc/2 

m 2fc+1 = M 2k + 2 J 2n+1 G ( 2n + 1, 2fe + 1) 

n=l 

where we have introduced a dimensionless parameter J7" = = ^tj^ and the 

?77q 1V1 

function G(l,h). For every coefficient rrih this function describes the numerical 
factor multiplying the power / of the parameter J . 

By substituting the expressions (4) in the expansion (2) of the Ernst's po- 
tential £ on the symmetry axis, and rearranging sums, it is possible to write this 
potential as a power series of the parameter J: 

M M 2 A 

Z(p=0, Z ) = -+J— + Y,J a ®a , (5) 

ce=2 
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where functions $ Q are defined below 

oo 

$ 2n = G(2n, 2k)X 2k+1 (6a) 



k=2n 
oo 

$2n+i = Yl °( 2n + 1 > 2k + l)A 2fc+2 , (66) 

k=2n 

M 

with the notation A = — . Let us note that since the parameter J is imaginary, 

z 

the functions $ a with an odd index turns the series (5) into an imaginary function. 

As shown in (5), the Ernst's potential £ can be expressed on the symmetry 
axis by a double series; one of them in terms of the parameter J and the other 
being a power series of the inverse coordinate z. Nevertheless, as we will see, the 
sum of the series <E> Q can be obtained, at least up to first orders. In order to do 
that it is necessary to obtain the analytic expressions which describe the double 
index functions G(a, h). If the first of those indexes is fixed, that is to say, if we 
consider a certain value for the power of the parameter J7", we have tried to adjust 
the resulting series of the corresponding terms arising from every coefficient m n . 

For example, it is very easy to check that factors appearing with powers two 
and three in the parameter J verify respectively the following expressions 

15 

G(2,2A:) = (2jfc + 3)(2fc + l)(2jfc-l) 

(7) 

G(S 2k + 1) = 15(10fc-17) 
K ' ' (2fc + 5)(2£; + 3)(2fc + l)(2fc-l) 

Now, it is quite simple to obtain the sum of the series $ Q by rewriting the functions 
G(a, h) as a sum of irreducible fractions. Particularly, for G{2, h) we have 

G(2,2j + 1) = 

,( 2 ) (8) 



1 ^ 

G(2,2j) = ^ 



1 2 1 

+ 
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o 2i + ^ 



2j + 3 2j + l 2j-l 
Taking this expression into (6a) the function $ 2 gives 

oo 2 (2) 
fc=2 j=0 J 

Rearranging sums and making use of Lemma 3 of Appendix, we obtain the fol- 
lowing finite sum: 

2 3-j 

$2 = ^^ 2) A 4 ^C 2( 3- j))2fc Q 2fc (l/A) , (10) 

j=0 k=0 
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where the coefficients Cih are defined in Appendix and the functions Qh{x) are 
special Legendre's functions of second kind. Let us note that the previous expres- 
sion can be written as follows 

3 min(3 — k,2) 

$2 = A 4 ^g 2fc (l/A) J2 9j 2)C 2(3-j),2k ■ (11) 
fc=0 j=0 



3. SEQUENCE OF EXACT SOLUTIONS 

The expression of Ernst's potential on the symmetry axis can be used as 
boundary condition to obtain solutions of the Ernst's equation. For example, Sib- 
gatullin's method [2] simplifies this problem by solving a linear system of integral 
equations. So, the Ernst's potential E results from the following expression 



E=- [ ^mda , (12) 

J-l 



1 mfrMr) 



where r is a complex variable defined from the cilindrical Weyl's coordinates 
t = z + ipa, and a G [—1, 1] is an arbitrary integration variable. Function e(z) 
represents the value of the Ernst's potential E on the symmetry axis, that means, 
e(z) = E(p = 0,2:). At length, the function p(a) must be a solution verifying the 
following integral equations system 

p C h ±^M da = (13a) 
J -i (r — 77) vl — a 2 

1 f 1 fa 



K J -1 y/i — a 



da = l , (136) 



being rj a complex variable defined as i] = z + ips, with q e [—1, 1], and where the 
symbol p stands for the taking of the principal part of the integral. On the other 
hand, the function h(r, rj) is defined as is shown below 

h{r,ri) = e(T) + e(ri) , (14) 

where the function e(rj) is obtained from e(rj) by conjugating first the variable, 
77 — > rj* , and then the function, i.e., e(rj) = e*(rj*). 

Obviously, the general solution of the equations in (13) is not evident. Never- 
theless, rather compact expressions have been obtained [9] for the Ernst's potential 
when the boundary condition e(z) is a rational function, i.e, 

E(p = 0,z) = e(z) = ^ , (15) 
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where P(z) and Q(z) are polynomials in the variable z, which, taking into account 
that the Ernst's potential must tend to 1 in the neighbourhood of infinity, should 
be as follows 



N 



P(z) = Z N + J2a k z 



N-k 



k=l 
N 



(16) 



N-k 



k=l 



Since we know the structure of the potential £ on the symmetry axis in terms 
of coefficients m n , inmediately a question arises: wether there exists a relation 
between these coefficients and those of the polynomials in (16) . The above question 
has been answered in [11] and in what follows we give the resulting expressions 
for P(z) and Q(z) in terms of m n : 



P(Z) = (I/Tv) -1 



N-i _ sr^ N ~ 2 



y N-2-n 



Q(z) = (L N ) 



-1 



z — mo 
1 



* + L„=0 m nZ 



y N-l-n 



Z + En=0 m n Z 



z + m 
1 



y N-2-n 



m N 



m N -i 



mi 
m 



m N 



m N -x 



mi 
m 



m 2 N-i 
m 2 N-2 



m N 
m N -i 



m 2 N-i 



m 2 N-2 



m N 
m N -i 



(17) 



(18) 



where the following determinant Ln has been defined: 



'N 



m N -! 
m N - 2 

m 



m N 
mjv-i 

mi 



m 2N - 2 
m 2N - 3 

m N -i 



(19) 



Hence, from 2N coefficients m^, it is possible to build on the symmetry 
axis the Ernst's potential E of a vacuum solution, which is a ratio of order iV 
polynomials. 

Let us consider now the following question: is potential E of the Monopole- 
Dynamic dipole solution a polynomial ratio on the symmetry axis?. On handling 
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the coefficients m n in (3) corresponding to such solution, the determinants L n 
(19) seem to be unlikely to be zero originating from some order iV onwards. Since 
the coefficients m n are not available for every order, we can only assert that the 
behaviour of the Ernst's potential of such solution does not correspond to a ratio 
of polynomials of order N < 10. 

In spite of last statement, it is possible to construct a set of rational type 
potentials E on the symmetry axis involving these coefficients m n . Thus, by 
Sibgatullin's method one can obtain a sequence of exact solutions which, as will 
be shown, approach the Monopole-Dynamic dipole solution. 

The coefficients m n (3) have been obtained on condition that the multipole 
moments higher than the dipole are zero. Hence, on fixing N coefficients m n , 
we will get an Ernst potential which describes a solution whose N — 2 multipole 
moments higher than Angular momentum are zero. At the same time, multipole 
moments of higher order, although different from zero, are determined by just 
those N coefficients mu- 

In order to perform the sequence of exact solutions, let us proceed to con- 
sider the Ernst's potential on the symmetry axis as a ratio of polynomials whose 
order N will be progressively increased. Hence, at each stage, we will be fixing 
an increasingly bigger 2A number of multipole moments for the corresponding 
solution. 



A) ORDER N = 1 

Let us suposse first that Ernst's potential on the symmetry axis is a ratio of 
polynomials of order N = 1, i.e., 

(1) z + a 1 _ PW(Z) 

6 {Z) ~ z + b 1 ~ Q(D(z) • {M) 
In order to calculate coefficients cii and bi we handle the two first coefficients nik 
in (3), 

m = M , rrn=iJ . (21) 

We obtain polynomials P^ 1 \z) and Q^(z) in terms of these coefficients by using 
expressions (17) and (18) respectively, and as a result the Ernst's potential is 
written on the symmetry axis in the following way 

e [Z) - z + M-iJ/M ■ [ZZ) 
The previous expression is exactly the corresponding potential of Kerr's metric 
with parameters M and a = J/M. By using Sibgatullin's method and the expres- 
sions in [11], it is possible to obtain the Ernst's potential for every range of Weyl's 
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coordinates {p, z}, as follows 



E Kerr = a ( r + + r ~) + ia ( r + ~ T '-) ~ 2aM ^ 

a(r+ + r_) +ia(r + -r_) +2aM ' 1 ' 

being a the positive root of polynomial P(z)Q(z) + P{z)Q(z) (numerator of the 
function h(z,z) (14)), i.e., 

a± = ±^M 2 - a 2 (24) 



and where r±= a/ p 2 + (z — a±) 2 

Now, let us calculate the coefficients rrik higher than those proposed in (21). 
To do that, we can use some expressions in [11] which relate such coefficients to 
those of the polynomials P^^z) and Q^'(z), and gives 

m k = M{ia) k = M k+1 J k . (25) 

It must be remembered that one property of Kerr's metric turns out to be 
the identity between its multipole moments and the coefficients rrik entering the 
expansion of the potential £ on the symmetry axis [1], and leads to, 

M = M , M 1 = JM 2 , M 2 = J 2 M 3 , M 3 = J S M 4 , 

M 4 = J 4 M 5 , M 5 = J 5 M 6 , M 6 = J 6 M 7 . 

(26) 

Obviously, the coefficients m n higher than mi do not equal the corresponding 
coefficients of the M — J solution. That is a good reason to step forward. 

B) ORDER N = 2 

Let us consider the Ernst's potential on the symmetry axis as a ratio of 
polinonials of order iV = 2, i.e., 

(2) _ z 2 + g lZ + a 2 N^(z) 

[Z) ~ z 2 + b 1 z + b 2 ~ DW(z) ' { <] 

We introduce four coefficients m^, according to expressions (3) with the following 
values 

mo ee M , mi ee ij , m 2 = TO3 = . (28) 

If we choose the coefficients rrik in (28), the solution we generate will have the 
quadrupole moment and the octupole moment equal to zero since m 2 and ms are 
just equal to these multipole moments respectively, . 
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The calculations of the polynomials P^ 2 \z) and Q( 2 \z) leads to the following 
result 

MZ ~' J (29) 



: (2) M = - 



z 2 +Mz + iJ 

In order to construct the Ernst's potential according to Sibgatullin's method 
it is necessary to obtain the roots of the function h(z, z) = P(z)Q(z) + P(z)Q(z). 
For this case, two roots of that function turns out to be real numbers whereas the 
other two are imaginary conjugated numbers: 



M 



(30) 



^ v - v - -~ 2 
From these roots it is possible [11] to write out the Ernst's potential E as follows 

A + r 



E^(p,z) 



A-r 



(31) 



A = -r{ 



+ 



p — 1 



P 



+ 



p + l 



V 



1 

+ -r 1 



P~ 1 \ + i 

— r 2 + 
P 



P+l 



P 



M 

T 



(p-p) 



p — 1\ ( p + 1 



v 



p 



r + _ 



l p 2 p 2 - 1 . _ + 

4 p 2 p 2 11 22 
p + 1 \ f p — 1 



p 



P 



M 

+ yb + W 



p — 1 \ f p — 1 



P 



P 



p + l\ (p+l 



p 



P 



where the following notation have been used 



p = +Vl + 2 J 



^ = +Jp 2 + (z-af) 2 



(32) 



(33) 



Oddly enough, for this case the structure of the potential £ on symmetry axis 
is defined by the following coefficients m n 



mo = M , mi = iJ , rrik = , V/c > 2 . 



(34) 
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and so, reproduce only up to order 4 the coefficients that characterize the M — J 
solution. Multipole moments of this solution turn out to be 



M = M M6 = ~h j2M? 



M 1 = JM 2 M = ™_j* M i 

429 

M 2 = , , 

- - - ,n / J- 



M3 = «> = M " (w^ - 3^3 ^) (35) 

10 V 4199 12597 

It should be noticed that quadrupole moment and octupole moment are zero 
by contruction. The first multipole moment different from zero is M4, which turns 
out to be proportional to J 2 , one order higher than the angular momentum in the 
parameter J . It is noteworthy that all higher massive moments are proportional 
to J 2 and the dynamic moments turn out to be of order 3 in that parameter, and 
therefore, just the same order than quadrupole moment and octupole moment in 
the previous case N = 1. We will discuss this issue later. 

Another interesting property of this solution is its equatorial symmetry which 
can be inferred from the fact that its odd multipole moments are imaginary quan- 
tities (and so, according to notation FHP [1], represent dynamic moments ) 
while even multipole moments are real quantities (masive moments). It is an in- 
trinsec characteristic of the procedure used to construct solutions from rational 
type Ernst's potentials on the symmetry axis. In fact, it can be proved that if 
coefficients m n introduced are alternatively real and imaginary quantities, then 
the resulting Ernst's potential has equatorial symmetry. That occurs because the 
coefficients and bk (16) fulfill the next relation at = (— l) fc b* k . Hence, in order 
to have equatorial symmetry an axisymmetric stationary and asymptoticaly flat 
vacuum solution must meet the following necessary sufficient condition [12], [13]: 

e+(z) e* + (-z) = 1 , (36) 

where e+(z) denotes the Ernst's potential on the positive region of the symmetry 
axis and symbol * denotes complex conjugation. 

C) ORDER N = 3 
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In this case we introduce six coefficients m^ according to the ones possesing 
the M — J solution (3), which means, 



mo = M , mi = iJ , m,2 = ms = , 

1 2 Is ( 3? ) 

m 4 = -MJ 2 , m 5 = -i^-J ■ 

This choice ensures that the multipole moments lesser than M 6 are zero, except 
for mass and angular momentum. 

Then, the Ernst's potential on the symmetry axis can be written as follows 

e^(z)-^l (38) 

being 

P^\z) =z 3 + Mz 2 (-1 + J/3) + zM 2 (-l/7 - 4J/3) + 

+ M s (l/7 + J/7-J 2 /3) / 

(39) 

Q^\z) =z 3 + Mz 2 (l + J/3) + zM 2 (-l/7 + 4J/3) + 
+ M 3 (-l/7-J/7-J 2 /3) . 

The roots of the function h(z, z) (14) corresponding to this potential are 



af =±M X /A-+A+ + ^ + ^ 



af = ± M^- l -(A- + A+) + 1 + ^ + i\V3(A+ - A~) , (40) 



af =±M^-^(A-+A+) + ^ + ^- i\V3(A+ - A~) 
where the following notation has been used 



A ± = (a± V&) (1 



/3) 



8 106 2 131 4 1 6 
a= 343 "441 J + 3402 J + 19683 J ( 41 ) 
5 = _ 72 , 21284 T 4 , 71546 6 5309 8 1 10 

50421 583443 6751269 4960116 177147 

If |J7"| < 1, which is realistic enough, then b > and a ± Vb > 0, and so, there 
are two real roots and two pairs of complex conjugated roots. For the sake of 
concision we will dispense with the expression resulting for the Ernst's potential. 
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The coefficients m k of this exact solution can be calculated by using expres- 
sions in [11] to obtain the following ten coefficients 



mo =M , mi — JM , mi — — , 
m4 =-iM 5 J 2 , m 5 = i«V 

™,=M«(^ + ^) (42) 

According to these expressions the multipole moments are 

M = M , Mi = M 2 J , M 2 = , M 3 = , M 4 = 

M 5 = , M 6 = -M 7 (-ii7^ 2 + el^ 4 ) > M ? = M 8 (--^ 3 + ^J 5 ) 

A « - M9 <-ai J2 + ^ ~ Hot J6 » 

,.m/ 26912 o 9158 ^ 1 ,, oN 

M 9 = - M 10 ( J 3 H J 5 J 7 ) (43) 

v 2501499 3216213 1701 ' 

wii/ 13392 ^2 55500 , 130 . 1 

Mio = - M n ( J 2 H ,7 4 J 6 H J 8 ) 

u v 10081799 15842827 1281987 5103 ' 

We can see that, by construction, this solution obviously posseses a higher number 
of null multipole moments than the previous solution. Besides, the first multipole 
moment different from zero, i.e., Mq, turns out to be proportional to J 2 again, 
that is to say, one order less than the angular momentum. Nevertheless, it must 
be pointed out that its magnitude is not necessarily smaller than the first moment 
different from zero (M 4 ) in the previous case. In fact, first multipole moment 
different from zero for each case is always proportional to J 2 , and so, we can achive 
striking solutions in this process which, in comparisson with previous solutions in 
the sequence, possess moments of a higher multipole order and higher magnitude 
at the same time. 

In order to illustrate the behaviour of the multipole moments in this sequence 
of exact solutions, we write out the moment of a certain multipole order for each 
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solution. By way of example, we will compare the moment M 6 of each solution 
with just the same moment for the case N = 1, i.e., Kerr's metric 



= M 7 J 6 = M* err 

M (2) - L/Lf 7 T 2 - 1 -M Kerr — f A A\ 

M 6 - J - 33 M 6 ^ (44) 

M< 3 > = --M 7 J 2 - -M 7 J 4 = - — M* e ^ 

6 49 63 49 6 J 4 63 6 J 2 

That is to say, the higher approximation degree in the series of solutions, the higher 
magnitude of the multipole moment. According to (44), the sequence of solutions 
should have a good behaviour, i.e., progressive diminution of the magnitude of any 
moment, if the parameter J were larger than 1 (absolute value), i.e., J > M 2 , 
which is not an expected condition for any realistic object. 

For these reasons, and others which will be disscussed in next section we 
will introduce a different approach to the Monopole-Dynamic Dipole stationary 
solution. 



4. STATIONARY APROXIMATE M-J SOLUTION 

According to the previous section, the construction of an exact stationary 
and axisimetric solution by Sibgatullin's method requires the structure of Ernst's 
potential on the symmetry polynomial ratio. Nevertheless, we have shown 

that the structure of the Ernst's potential on the symmetry axis corresponding 
to a solution of type M — J is, in some way, a ratio of series, which means that 
it cannot be expresed as a polynomials ratio. Hence, although the potential of 
such a solution were obtained on the symmetry axis, we cannot apply Sibgatullin's 
method. Besides, the sequence of exact solutions previously proposed approches 
the M — J solution in a rather unexpected way, since the magnitude of multipole 
moments does not decrease while the approximation degree raises. 

Therefore, in this section we will proceed to approach the M — J solution in 
a different manner. We give up looking for exact solutions and propose instead a 
sequence of approximate solutions as partial sums of power series on the parameter 
J. 

The expressions in (3) of the coefficients m n corresponding to the solution 
M — J lead to an Ernst's potential £ on the symmetry axis as a power series in the 
parameter J . According to this result we will look for solutions in that way. Let 
us consider the Ernst's equation for the potential £ and let us assume a solution 
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of the form 



oo 



^ = ^ + J2^aJ a , (45) 



Q = l 



where £ represents the the Ernst's potential corresponding to Schwarzschild's 
solution. Imposing this series to verify Ernst equation at each order leads to the 
following equations concerning the functions £ a 



(£ 2 - l)A£ 2 a+i - 4£ V£oV£ 2 a+i + 2£ 2a +i(V£ ) 2 = H 2a+1 
{il - 1) A£ 2a - 4^ V6V6a + 26 a (v^) 2 ^-^ = ^2 



£o 2 -l 



where first equation refers to odd orders and second one to even orders (a = 
1,2,...); the second members of those equations are given by 

H a = (-l)M26V^Va-^A^] , a>0 . (47) 

i-\-j-\-k=ot 

i,j,k<a 

i.e., equation of order a depends on the previous orders. 

The previous equations can be simplified by redefining the functions £ a as 
follows 

which leads at each order a to the following equations 



Ca = ^ • (48) 



(^0 2 -l)ACa+^ 2 Ca = ^ I , <( (49) 



v = , a = even 
v = 2 , a = odd 



It is easy to solve this equation by writing it in prolate coordinates. As usual the 
general solution can be obtained by adding a particular solution of inhomogeneous 
equation to the general solution of whole equation. Moreover we impose a regular 
behaviour on the symmetry axis (y = ±1) at least like 1/x with respect to the 
variable z in the neighbourhood of infinity. We then obtain: 



x 

Qf=l 



1 OO 

(50) 

C^X,y) + J2 h nQ i n ) ^)Pn(y) 



X 2 -1 

£,a{x,y) = 



X 2 



n=0 
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where Qn{x) are associated Legendre's functions of second kind, the functions 
Ca{ x i v) are particular solutions of the inhomogeneous equations corresponding to 
each order a, and h% are arbitrary constants. 

To describe the Monopole-Dynamic dipole solution (M — J) from the general 
solution in (50), it is necessary to add as a boundary condition the behaviour 
of potential £ on the symmetry axis, which has been defined previously by the 
series in (5). Hence, we force now the function $ a appearing in (5) to agree with 
the corresponding restriction on the symmetry axis of the function £ a of general 
solution (50), which leads to determine the constants h%- 

Previously, the functions $ Q must be adapted to the structure of the general 

M 2 - z 2 

solution in (50), and so, we begin by taking a factor out of the expression 

z 

of the function $ Q , 



M 2 - z 2 

$2n = 



Z 2 



1 oo 

^ G(2n,2k)X 2k+1 



^ 2 ' /, =2 n 



(51) 



By carrying out an expansion on the parameter A results in 

^n = -^^f^Y,G(2n,2k)X 2k + 2 ^ , (52) 

i=0 k=2n 

that is to say, 



$2n = -^^^A^ +1 ^G(2n,2fc) . (53) 

j=2n k=2n 

Below we write these functions $2n in terms of the Legendre's functions of 
second kind, by using Lemma 4 of the Appendix 

T\/t2 2 00 00 

*2n = £ lf n) £(4z + l)L 2i>jy g 2i (l/A) , (54) 

j=2n i=j 

where the following notation has been used 



oo 



lf n) = E G(2n, 2k) . (55) 



k=2n 



With respect to the odd functions $2n+i we proceed in the same way to obtain 
mentioned factor and write, in this case, these functions in terms of associated 
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Legendre's functions of second kind Qff +1 (l/\), which can be obtained making 
use of Lemma 5 presented in the Appendix, and so, 

r 2 00 ,7.0 I 



z 2 — M' x 



2a+l 



v 4Z + 3 (2) r ^ 21 + 2n + 1 (2q+1) 

2^ r9/-U9V9/-H) y 2Hl^/AJ 2^ ( 2n+ l)M ^2Z,2n^„ 



^(2/ + 2) (2/ + 
with the following notation 

n 

/( 2a+1 ) = ^ G(2a + l,2fc+l) . 



(56) 



(57) 



fc=2a 



At this point we proceed to determine the constants h™ of the general solution 
(50) which correspond to the M — J solution . In that way, we choose particular 
solutions of inhomogeneous equations (49) as follows 



C2 P a=E( 4Z + 1 )^(^)^(y) 



2=0 



(58) 



C 2 P a + i =^{41 + 3)Q%> +1 (x)S 2l+1 (y) 



1=0 



S a {y) being polynomials in the angular variable. 

Comparing general solution (50), evaluated on the symmetry axis, with the 
previous expressions (54) and (56) gives 



h 2a 
n 2n+l 



^ = -(4n+l)S 2n (l) , n<2a 







^2n = -(4n + l) 



hl a n +1 = 



h\nX\ = -(4n + 3)S , 2n+i(l) , n<2a 



S 2 n(l) + J2 

L2n,2kl, 



(2a) 
k 



k=2a 



n > 2a 



h 2a+l 
a 2n+l 



= -(4n + 3) 



5 2 „+i(l) + 



1 



E2n + 21 + 1 r (2a+i) 
— . s.. L/9n..9.lJl 



(2n+l)(2n + 2)^ a (2/ + 1)!! 



n > 2a 
(59) 

Let us construct explicitely the first orders of the solution M — J. Obviously, 
the order zero contribution to the solution must be Schwarzschild's solution, since 
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taking the parameter J = leads us to consider the mass as the unique multipole 
moment. The Ernst's potential £ of Schwarzschild's solution (£ = 1/x) equals 
the structure described on the symmetry axis (5). 

A) FIRST ORDER. 

The first contribution on the parameter J should be a solution of the first 
equation in (46) turns out to be homogeneous at order one, and hence 



i 2 00 

Hi = -^-J2 h iQi 2) W p ^ ■ ( 6 °) 

Jb 

1=0 



This expression on the symmetry axis, gives 

< 

£ 1 (y = l)=2A^+2A 2 /ii + (A 2 -l) 2 ^/i, 1 Qf ) (l/A) . (61) 

1=2 



In addition, the first contribution to the M — J solution must agree with (5) on 
the symmetry axis and so, the only solution corresponds to the following choice 
of constants 

hh = , h\ = \ , hl=0, l>2 , (62) 

which means, 

U^y) = l 1 —^Q i " ) (x)P 1 (y) = ^ . (63) 

Zi Jb Jb 

It must be noted that this approximate solution is the same as the one arising 
from the expansion of Kerr's metric on parameter J up to the first order. 

B) SECOND ORDER. 

A particular solution of the inhomogeneous equation (46), corresponding to 
order two, results in 

= ( I " T>T^ ' (G4) 

and in view of Lemma 4 of the Appendix, this could be rewritten in terms of 
associated Legendre's functions in the following way 



C 2 nh = J2Q2n(x)(4n + l) 



n=0 



- -y 2 (l -L 2n ,o) 



(65) 
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By substituing this particular solution in the expressions (59) we determine the 
constants of the contribution of the second order to the solution, i.e., 



1 



^0=0 ' ^2 = 5 L 2,0 ~ 7T 



1 



h 2 2n ={An + l) 



1 n 

~ + L 2n fi - ^ L 2n,2kl k 



k=2 



(66) 



, n > 2 



(2) 

Now, taking the expression G(2, 2j) (8) into account we can write term I k ' as a 
sum of irreducible fractions in the following form 



f (2) = (fc + 3)(fc-l) 1 15 1 15 1 
k (2fc + 3)(2fc+l) 4 8 2/c+l 8 2k + 3 



(67) 



which implies that 1^ = —1 and l[ 2 ^ = 0. Therefore constants h 2n with n > 2 
can be written as follows 



h\ n = {An + 1) 



1 

-- - ^2L 2n ,2kI l 



(2) 



fc=0 



(68) 



Considering now expression (67) we have 



h\ n = {An + 1) 



1 15 \ ^ L 2n ,2k _ 15 ^ L 2rii2k ^ 1 \ ^ 

2 8 ^2k+l 8^2^ + 3 4^ 2n ' 2fc 

fc=0 k=0 k=0 



(69) 



and making use of Lemma 2 of the Appendix, and the ortonormality of the Leg- 
endre's polynomials, we have finally the following expressions 

h 2 = \ , h 2 = -5 , fc| n = -§(4n+l), n>2 . (70) 

Hence, the contribution of order 2 to the solution M — J can be written as follows 

1 ^ ' 1 " ^ \Qo{x)P {y) - 5Q 2 {x)P 2 {y)- 

(71) 



&(x,y) = 7; 



2x x 3 



or 



- J]-(4n + l)g 2n (x)P 2n ( 2 /)] , 

n=2 

an expression which can be simplfied by making use of Heine's identity [14]: 

00 

Y,(^ + l)Q 2n {x)P 2n {y) , (72) 
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x 

x 2 — y 2 



n=0 



which leads to the following expression 



, (73) 



that is to say 



5 x z 



4 x 2 



Qo(x)P (y) - Q 2 (x)P 2 (y) 



X 



x 2 — y 2 



. {Kerr 
+ ?2 



(74) 



where ^ err is the order 2 in the expansion of Kerr's metric for the parameter J . 

The multipole moments can be calculated by FHP algorithm from coefficients 
m n . For the solution up to order 1, i.e., £m_j = io + J £i 5 all coefficients m n with 
n > 2 are zero (i.e., of higher order than J\ So, its multipole moments equal 
those in the exact solution presented in the previous section whose potential e(z) 
was a ratio of polynomials of order 2, that is to say, 



M =M , Mi = JM' 



M 2 = , M 3 = , 



1 



M 4 =^J 2 M 5 
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M 7 = J^M 8 

7 429 



M 5 = -— J 3 M 6 



M s = -AT 



M 6 = - — J 2 M 7 
5 33 



1 



:J 2 + 



53 



M 9 =M 



10 



43 41 - 



2431 



17017 



143 3003 

, M 10 = M 11 



4199 



7 ^ 2 202 < 



12597 

r (2) 



(75) 

With respect to the M — J solution up to the second order, i.e., Cm-j = £o + 
i7 £i + J" 2 £ 2 7 we have the following multipole moments 

M = M , Mi = JM 2 , M 2 = , M 3 = , M 4 







4A/T9 



J 4 M 



M 9 = M 



10 



V17017 51051 



M10 = 



49873 
6789783 



(76) 

It can be seen from (75) and (76) that the higher the order of approximation 
to M — J solution, the higher (one order more) the order in the parameter J of 
its non vanishing multipole moments. 

The structure of coefficients rrik in terms of the parameter J shows which 
is the order n of such coefficients that possess a contribution of order a in J . 
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In fact, the first contribution to an even power in the parameter, arises from the 
coefficient of order 2a. If a is odd, such contribution arise from the coefficient of 
order 2a — 1. Since the relation between the multipole moment and the coefficient 
mfc of same order is linear we can conclude that the multipole moments of the 
solution which approche the M — J solution up to the order a in the parameter 
J have the following characteristics: 

1) If the order a is even, all its massive multipole moments up to M 2a +2 
(incl.) will be zero, the following ones being at least of order a + 2 in parameter 
J . With regard to its dynamic multipole moments they will be zero up to M 2a -i 
(incl.), and the following ones will be at least of order a + 1 in parameter J. 

2) If the order of approximation a is odd, then all massive multipole moments 
will be zero up to M 2a (incl.) and the ones higher than that will be at least of 
order a + 1, whereas the dynamic moments will be zero up to M 2ce+ i and the 
following moments will be at least of order a + 2. 

These results show, in the same way as the static M — Q solution [5] , that it 
is possible to understand the series in (45) in terms of the perturbations theory. 
Each partial sum of that series is a better approximation to the solution which 
only has mass and angular momentum, since the multipole moments higher than 
this are either zero or have an order in the parameter J higher than the one 
of the approximation. In addition, unlike the solutions in the previous section, 
the multipole moment of a certain order is progressively smaller as the order of 
approximation in the series (45) grows. 

APPENDIX 

In this Appendix we will enunciate a sequence of Lemmas about several prop- 
erties of the Legendre's polynomials and the associated Legendre's functions of 
second kind. Results of these Lemmas are probably well known, but proofs for 
them can be easily obtained considering some results of various Lemmas appearing 



Let us P 2n be a Legendre's polynomial of even order in an arbitrary variable 



in [5]. 



n 




(A.1) 



where its coefficients have the following expression 
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Let us consider the development of an arbitrary variable to an even power in terms 
of Legendre's polynomials in that variable, i.e., 

00 

C 2n = E C 2n,2A(C) , (A3) 
k=0 

where coefficients C 2n ^ 2 k can be obtained by integration from the following ex- 
pression, 

Ah _i_ 1 r 1 

C 2n ,2k = — ^- / P2k(()( 2n d( , (A4) 



2 



1 



and so, for several values of indexes k and n such a coefficients turns out to be 

2n! 

C 2n , 2k = { {4k + 1) (2n-2k)!!(2n + 2k + l)H ' * ~ U . ( A5) 

: k > n 



Lemma 1. The following orthogonality relation is satisfied: 

fc 

E L 2 k, 2 jC 2 j j2n = 8kn • (A6) 
j=0 

Corolary: It can be deduced evidently another orthogonality relation 

i 

E L 2 k,2nC2j,2k = 8j n ( A7) 



Lemma 2. For all pair of positive entire numbers n and k such that n < k, the 
following equality is verified 

y - L2k ^ - = . (as) 



i=o 



Lemma 3. For all pair of positive entire numbers a and j the following equality 
is verified: 



E 2n + 2l + 1 = ^ E ^ + 4a,2nQ2n(l/A) , A EE - . (A9) 
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2n + 2j + _ 

n=2a J n=0 



Lemma 4. For all positive entire number n following is verified 

^ oo 

k=n 
^ OO 

= (4H3)L 2H1 , 2)1 _ 1 Q 2H1 (x) . (A.lOb) 



x^ n ~ 



k=n—l 



Lemma 5. For all positive entire number n following is verified 

E (^C.!*^^ • (AH) 

fc=n — 1 
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